Antibunching correlations in a strongly coupled exciton - photonic crystal cavity 
system: Role of off-resonant coupling to multiple excitons 
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We employ a master equation approach to study the second-order quantum autocorrelation func- 
tions for up to two independent quantum dot excitons, coupled to an off-resonant cavity in a photonic 
crystal - single quantum dot system. For a single coupled off-resonant exciton, we observe novel 
oscillatory behavior in the early-time dynamics of the cavity autocorrelation function, which leads to 
decreased antibunching relative to the exciton mode. With a second coupled exciton in the system, 
we find that the magnitude and the lifetime of these oscillations greatly increases, since the cavity 
is then able to exchange photons with multiple excitonic resonances. We unambiguously show that 
this spoils the antibunching characteristics of the cavity quasi-mode, while the autocorrelation of 
the first exciton is unaffected. We also examine the effects of detector time resolution and make a 
direct connection to a series of recent experiments. 



Introduction.- Two level behavior is ubiquitous in 
quantum physics, and when a two- level system is cou- 
pled with a suitable optical cavity system, the regime 
of cavity - quantum electrodynamics (cavity-QED) can 
be realized. Cavity-QED is a large and fascinating field, 
and several well known effects have been observed and 
exploited for a number of years, e.g., the Purcell ef- 
fect [l[ and vacuum Rabi splitting 0. Cavity-QED sys- 
tems that emit deterministic single photons are also in- 
teresting from an applications viewpoint, opening doors 
to practical ways of doing quantum information process- 
ing 0, 0|. Two-level system behavior can occur in a 
number of material structures, including atoms 0, @], 
ions 0], molecules [1, [p , and color centres [Icj. Semi- 
conductor QDs fnl.ll2l| have also been proposed as an 
"artificial atom," with strong excitonic levels that mimic 
two level behavior and are particularly attractive since 
they can function as a scalable and compact system for 
emitting single photons. Direct coupling between QD 
excitons and semiconductor cavities has recently been 
demonstrated in a variety of s yst ems such as micropil- 
lars (l3l.fT3.fl5 l, m icrodisks [THH?]], and photonic crystal 
(PC) cavities [il [H. 

An important difference between the modern QD semi- 
conductor devices, and earlier atomic systems, is that the 
confined photon cavity environment is considerably more 
complicated than a couple of mirrors, frequently exploit- 
ing photonic bandgap physics and high-indcx-contrast 
modes. In addition, QDs exhibit rich excitonic struc- 
ture, hence it is not known to what extent the two level 
approximation applies, though it is expected to be suit- 
able for narrowband cavity systems with a well defined 
cavity resonance that couples to the target exciton mode. 
Thus, it can be anticipated that, in addition to the fa- 
miliar cavity-QED phenomena well known to atoms, fur- 
ther complexity unique to these semiconductor cavity - 
QD systems can been observed. A recent example of 
off-resonant coupling between a cavity and a single QD 
resulted in the observation of a pronounced cavity mode 
emission [l^, [2(| [H| ■ First considered a surprising result, 
it is now known that this is likely due to the non-trivial 



coupling between the leaky cavity mode and the exci- 
ton mode [13, [H, [24]. J25I [26[ , which is exasperated in a 
planar PC system [261 ] - Even more surprising has been 
the various reports of poor antibunching - a measure of 
the quantum nature of the emitter - even in the strong 
coupling regime for an exciton detuned from the cavity 
resonance [191 ]; and further experimental studies of the 
quantum correlation function in the off-resonant regime 
have provided additional perspectives on these results. 
Specifically, it has been observed that the cavity autocor- 
relation function shows significantly worse antibunching 
than the exciton autocorrelation function [2(| , and to the 
best of our knowledge this perplexing behavior has not 
been yet explained. Yet an understanding of this mystery 
is a crucial step towards the understanding and design of 
QD-based single photon sources and would provide fresh 
modeling insight into the nanoscale light-matter interac- 
tions at a very basic level. 

In this Letter, we introduce a quantum optics formal- 
ism that allows one to calculate both the exciton and the 
cavity mode quantum autocorrelation functions for a sin- 
gle QD - PC cavity system. We provide a framework in 
which to understand the above experimental results, by 
emphasizing the consequences of the quasi-mode nature 
of the cavity, and the significant effect of additional ex- 
citons on the quantum autocorrelation functions of the 
system. We show the general need to account for cou- 
pling to additional excitons in the system that play a 
qualitatively important role, even if they are far off reso- 
nance. We investigate the nature of off-resonant coupling 
and pure dephasing in detail, and make a direct and suc- 
cessful connection to recent experimental results. 

Theory and Model.- Our model system consists of a 
single QD embedded in a planar PC cavity, but is also 
applicable to other semiconductor cavity systems where 
the leaky cavity mode dominates the emission, e.g., the 
micropillar cavities. Due to the leaky nature of the cavity 
mode, both the exciton and the cavity are able to emit 
photons out of the plane (vertically). In Fig. QJa), we 
show the relevant energy level diagram, that consists of 
up to two excitons (both originating from a single QD and 
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assumed to be independent - though they are indirectly 
coupled through the cavity mode) and the resonant mode 
of the PC cavity. The uppermost exciton energy level is 
utilized for exciting the system with a far off-resonant 
pump pulse. We adopt a master equation approach (27j 
to calculate the second-order quantum correlation func- 
tions. The master equation for the leaky cavity and the 
radiative decay modes, in the interaction picture, is 

P =-tt[Hi,p] - n[a\a c p + pa [ c a c - 2a c pa].) 

- [li{auf> + P&U - pa + )-j' l (a- li pcro + aooP^ii)] 

i=l,2 

- ^ liuipuuP + pbuu — 2a m pa ui ), (1) 

i=l,2 

with the interaction Hamiltonian 

Hi = ih J2 9i(*i&t ~ + (<4 - °Z) , (2) 

i=l,2 

where Q(t) is a classical excitation pulse that excites 
the upper lying exciton level, \x u ), 2-f iu is the fast non- 
radiative decay rate, a a = &f<j~, and p is the density 
matrix operator. The boson operators a c and a\ are the 
photon creation and annihilation operators for the cavity 
mode, and af and a~ are the Pauli raising and lowering 
operators, respectively. The parameters of the system 
include g - the coupling strength between the cavity and 
exciton; and 2ji, 2k and 2jfi - the radiative, cavity and 
pure dephasing decay rates, respectively. 

Importantly, both the one and two time equations of 
motion can be calculated from the master equation [27| . 
and utilizing the quantum regression theorem , we ob- 
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FIG. 1: (Color online) System under investigation, (a) En- 
ergy level diagram for theoretical model, (b) Photonic crystal 
slab with embedded quantum dot, showing a background ra- 
diative decay process and a leaky cavity decay process, (c) 
Cavity mode emission spectra for one (dashed red) and two 
(solid blue) excitons. 



tain the second-order quantum correlation functions 

G£U*,* + t) cx (vt(tK(t + r)K(t + T)K(t)) (3) 
G%(t, t + t) cx (al(t)al(t + r)a c (t + r)a c (t)) , 

for the exciton and the cavity mode, respectively. The 
time integrated (? ( - 2 ' ) ( T ) i s then obtained by integrating 
over time t for the center pulse, and normalizing to the 
peak of one of the sideband pulses [28| . 

Calculated Spectra and Antibunching Correlations. - 
We optically excite the system up to a higher-lying 
level (x u ) from the ground state via Gaussian excitation 
pulses, ri(t) , with full width at half-max of 10 ps spaced 
1800 ps apart; the pulse amplitudes correspond closely to 
7r-pulses to optimally feed the high-lying level. The sys- 
tem is then allowed to quickly relax from |a; u ) to \xi) via 
a fast non-radiative decay rate of 2ji U = 0.4 meV. The 
cavity emission spectra calculated in the presence of one 
and two excitons can be found in the lower panel of Fig.Q] 
In the red (solid line) spectrum, a single exciton couples 
to the leaky cavity mode and the cavity emission spec- 
trum consists of both the bare exciton resonance as well 
as the bare cavity resonance. These one exciton results 
are consistent with the results of others [HI, [24|, [2f| H|| . 
Also, since the cavity mode emission characteristics en- 
tirely dominate the overall collected spectrum [2^|, we 
only show the cavity mode emission in the spectrum. As 
recognized, introducing a second exciton has a small but 
non-negligible effect on the spectrum. Analogously to 
the first exciton, the second one also "feeds" the cav- 
ity mode, and therefore has the effect of increasing the 
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FIG. 2: (Color online) (a) Autocorrelation function for the 
cavity (solid red) and exciton (dashed blue) for pure de- 
phasing rates of 7/ = 0, 10 and 20 fieV from bottom to top; 
7 = 1/ieV in all cases, (c) Zoom about r = region of (a), 
(b) Degree of antibunching - the ratio of the area under the 
r = peak compared to peaks at other times. Inset: Associ- 
ated spectrum - for an exciton on resonance with the cavity 
in the strong coupling regime. 
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weight under the peak at the bare cavity resonance as 
well as introducing a peak at its own bare resonance. 
Again, we emphasize that the cavity emission spectrum 
will dominate the emission characteristics, and contains 
resonances at both the bare cavity resonance as well as 
the resonances of the excitonic modes that feed it. 

In choosing model parameters, we aim to mirror those 
determined from related experiments [ijj [HI, and 
choose a cavity-exciton coupling of g = 0.075 meV, and 
half rates of radiative and cavity decay of 7 = 1 /zeV and 
k = 0.05 meV, respectively. To determine an appropriate 
pure dephasing rate, we calculate the rate from a rep- 
resentative amount of antibunching for an on-resonant 
exciton and varying 7/. We then compare our results 
to those typically seen in experiment [2(| • The sec- 
ond order correlation functions are shown in Figs. [2a) 
and (c). The theory utilized is amenable to projecting 
onto either the cavity quasi-mode or the excitonic mode 
separately - even though the system is on resonance - 
which is something that is not experimentally accessi- 
ble. However, since the cavity mode dominates the emis- 
sion, the on-resonance condition is heavily dominated by 
the cavity mode contribution. With the chosen parame- 
ters, we extract a representative pure dephasing rate of 
7/ = 2 fieV - which corresponds to an antibunching dip 
of 23% on resonance (see Fig. [5] (b)). 

We now introduce a detuning between the exciton and 
the cavity resonance A xc = lj c — lu Xi where oj c and u> x are 
the bare cavity and exciton resonances, respectively. We 
examine a system with a singly coupled, off-resonance 
exciton x — x\, as well as an upper loading exciton, x u . 
The calculated autocorrelation functions are shown in 
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FIG. 3: (Color online) (a),(c) Cavity (solid red curve) and 
exciton (dashed blue curve) autocorrelation functions with a 
zoom in on the region near r = in the lower panel and 
the long time behavior in the upper panel. The detuning 
from bottom to top are A xc — 0.4, 0.8 and 1.2 meV. (b) Anti- 
bunching dip for the cavity (solid red) and the exciton (dashed 
blue). Inset: cavity emission spectrum schematic. 



Fig. [3] for three different detunings. In general, we ob- 
serve an overall broadening of all peaks (as compared to 
the resonant case) and a reduction in antibunching. This 
reduction is magnified with increased detuning. How- 
ever, for all values of detuning seen in Fig. [2] (b), we 
(2) 

note that g C c (t) shows decreased antibunching relative 

(2) 

to Qxx (t) , which agrees with the experimental observa- 
tions reported by Press et al. [20| , For larger detuning, 
the size of the antibunching dip levels off to about 62% 
and 67% for the exciton and the cavity, respectively. We 
note that the exciton is expected to experience some an- 
tibunching even in the absence of the cavity (28|, subject 
to the amount of pure dephasing, which would explain 
the limiting case seen for large detuning. A closer exam- 
ination of the early time dynamics of the autocorrelation 

functions for a single exciton (see Fig. E](c)) reveals new 

(2) (2) 
oscillations in the g cc ■ (t) that are not seen for g x J(r). 

These oscillations are caused by the off-resonant coupling 
of the cavity to the exciton, and are manifestations of the 
quasi-mode nature of the cavity mode. The period of the 
oscillations corresponds to the detuning frequency of the 
exciton T osc oc 1 //detune, and therefore the frequency of 
the oscillations increases with increased detuning. 

We next turn to describing the effects of other nearby 
excitons on the ensuing quantum statistics. We include 
a second exciton that is essentially identical to the ini- 
tial one (same 7, 7/, g) with the exception of its detuning 

with respect to the cavity mode, A X2C . We show a rep- 

(2) 

resentative plot of gcc (t) in the presence of one versus 
two excitons coupled to the cavity in Fig. 0] (a). For the 
case shown, the second exciton is doubly detuned from 
the cavity relative to the first exciton. In fact, in exper- 
iments, the coupling strength of the second exciton may 
be significantly larger than the first one (eg. Ref. fl9l |) 
which will only exasperate the results that follow. As can 

be seen in Fig.Q](a), the presence of the additional exci- 

(2) 

ton has profound influence on g cc (t) . These extraneous 
effects consist of magnifying all the features that were 
previously unique to g\c (t). In particular, we find that 
both the magnitude, as well as the lifetime of oscillations 
in the early time dynamics are significantly increased - 
thereby magnifying the differences in the antibunching of 
the cavity versus the exciton mode. 

In Fig.[H(b), we examine the effects of varying the de- 
tuning of the second exciton and summarize the results 
on antibunching. In this figure, we maintain an equal 
coupling strength of g — 0.075 meV for both excitons, 
and note that the oscillator strength of the second exci- 
ton line in the cavity emission spectrum decreases with 
increased detuning. We find that the antibunching of the 
stationary exciton (x\) is unaffected by the detuning of 
the other exciton - which is expected since we assumed no 
direct coupling between the excitons. For the second ex- 
citon (2:2), we find that increasing the detuning from the 
cavity resonance spoils the antibunching of the exciton, 
though it still retains a significantly stronger antibunch- 
ing dip than the cavity mode. For the cavity mode, we 
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FIG. 4: (Color online) (a) Cavity autocorrelation functions 
for one (smaller oscillations) and two (larger oscillations) ex- 
citons. The detunings for the two excitons are A xc — 0.8 meV 
and A xc = 1.6 meV. Inset: One (red) and two (blue) exci- 
ton spectra showing cavity mode emissions, (b) Antibunch- 
ing for the cavity mode (solid red), first excitonic mode 
(dashed blue), and second excitonic mode (dot dashed ma- 
genta). The detuning of the first exciton is constant at 
A X1C — 0.8 meV, while A X2C varies. Inset: Cavity emission 
spectrum schematic, (c) The effect of finite detector time 
resolution on g x 2 J(r) (top) and g x 2 x (r) (bottom) shown by a 
dashed black line. These were calculated for a single exciton 
that is detuned by A x±c — 0.8 meV from the cavity. 



find that in the presence of nearby additional excitons, 
the antibunching is largely destroyed, particularly rela- 
tive to that of the excitonic modes. However, as the 
second exciton is moved further off resonance the anti- 
bunching of the cavity mode improves, and for large A X2C 



we begin to recover the single exciton case. In cases where 
the second excitonic line is more intense than the initial 
one or with more excitons excited (e.g., [19(), we find 
that the substantial difference between the antibunching 
of the excitons and the cavity mode increases even fur- 
ther, again consistent with experimental observations. 

Finally, in order to connect more closely with experi- 
mental detectors, we include the effects of finite detector 
time resolution, by modeling the detector time response 
as a Gaussian function with a full-width at half-max of 
50 ps and an area of one. We apply this averaging to 



the output of the g^) (r) and gti (r) calculations and ob- 
tain the results seen in Fig. 4(c). The primary effect of 
finite detector response is an averaging out of any be- 
havior that occurs at time scales shorter than the 50 ps 
response time, such as the quick drop to zero seen on 
both the cavity and exciton autocorrelation for t — > 0, 

(2) 

and the oscillatory behavior observed in g cc (r). How- 
ever, this averaging preserves the relative areas under 
the peaks in the autocorrelations, and therefore has no 
effect on the size of the antibunching dips. 

Conclusion.- By utilizing a master equation model to 
describe the coupling between a leaky cavity mode and 
one or more QD excitons, we find significant differences 
between the second-order quantum autocorrelation func- 
tions of the exciton modes and cavity modes. These dif- 
ferences are manifestations of the quasi-mode nature of 
the cavity and include oscillations in the early time dy- 
namics of g cc ■ (t), as well as the reduced antibunching 
of the cavity mode relative to the exciton. Furthermore, 
we find that the presence of additional excitons magnifies 
this difference and largely destroys the antibunching of 
the cavity mode, which qualitatively explains recent, and 
hitherto unexplained, experimental observations [19|, l20j. 
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